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Perturbative coefficients grow factorially with the order and one needs a pre- 
scription to truncate the series in order to obtain a finite result. A common 
prescription consists in dropping the smallest contribution at a given coupling 
and all the higher orders terms. We discuss the error associated with this pro- 
cedure. We advocate a more systematic approach which consists in controlling 
the large field configurations in the functional integral. We summarize our best 
understanding of these issues for lattice QCD in the quenched approximation 
and their connection with convergence problems found in the continuum. 
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1. Introduction 

Perturbation theory has played an essential role in developing and establish- 
ing the standard model of electroweak and strong interactions. The renor- 
malizability of the theory guarantees that we can calculate the radiative 
corrections at any order in perturbation theory. On the other hand, a gen- 
eralization of Dyson's argument 1 suggests that the perturbative series are 
divergent and one needs to truncate the series. In absence of a definite 
prescription to deal with this problem, one usually relies on the "rule of 
thumb" which consists in dropping the smallest contribution at a given 
coupling and all the higher order terms. Clearly, this procedure has a lim- 
ited accuracy and it is not always obvious how to estimate the error or to 
decide if one needs to calculate one extra order. 

The problem is particularly acute for QCD corrections because they are 
large even at low order. As emphasized by Z. Bern's talk 2 , NLO corrections 
are important for multijet processes to be studied by the LHC. Another 
example 3 is the hadronic width of the Z° where the term of order a 3 is 
more than 60 percent of the term of order a 2 s and contributes to one part 
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in 1,000 of the total width (a typical experimental error for individual LEP 
experiments). It is not clear that the next term would improve the accuracy 
of the calculation. As these calculations can be extremely time consuming, 
it is necessary to address the lack of convergence of perturbative series in 
a more systematic way. 

2. The Rule of Thumb 

Consider a generic asymptotic series in a coupling A with coefficients grow- 
ing like CiCf T(fc + C 3 ). If we assume that for sufficiently small A the error 
made by truncating the series at order k is given by the order k + 1 contri- 
bution, it is possible to show that the error is minimized by truncating the 
series at order (A^l)^ 1 — C3 — (1/2). The rule of thumb leads then to an 
error which is approximately 

V27Ci(AC 2 ) 1/2 - C3 e-^ (!) 

This function is an approximate envelope for the accuracy curves at suc- 
cessive orders as illustrated in Fig. 1. 

Anh. Osc. 
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Fig. 1. Absolute value of the difference between the series and the numerical value for 
order 1 to 15 (in a Logio scale) for the anharmonic oscillator as a function of 1/A. As 
the order increases, the curves get darker. The dash curve is (log 10 (\/T2/w)c~ 3* ) 

This type estimate is not always correct. For instance, for the ground 
state of the double-well potential, the instanton effect is much larger 
than this estimate. Often, neither the asymptotic behavior nor the non- 
perturbative effects are known. Consequently, it is worth trying to figure 
out a general method to handle the problem. 
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3. Quantum Field Theory with a Large Field Cutoff 

As can be seen in several examples where a path integral formulation is 
available, the factorial growth of perturbative series is related to configura- 
tions with arbitrary large fields. Large order coefficients are dominated by 
large field configurations for which the expansion of the exponential of the 
perturbation at that order is not a good approximation. In two non-trivial 
examples 4 where the modified coefficients can be calculated numerically, it 
was shown that a large field cutoff drastically affects the asymptotic behav- 
ior of the perturbative series. At this point we are only able to do numerical 
calculations for specific types of field cutoffs. Namely, we remove the tails 
of integration in configuration space. If cf) x denotes a scalar field at sites x 
of a lattice, we impose the condition \(j) x \ < cf) max at every site. This pro- 
cedure is quite convenient for a Monte Carlo calculation, but not to write 
modified Fcynman rules because the additional condition is non-local in 
momentum space. At this point, a large field cutoff procedure that leads to 
simple Feynman rules remains to be found. 

Quantum field theory with a large field cutoff is a subject in infancy. 
Numerical studies 5 ' 6 show that, after appropriate rescalings, the transition 
between the large field cutoff regime and the small field cutoff regime can 
be described in good approximation by a universal function. This hints at 
a renormalization group explanation. It might be possible to construct an 
effective action, with couplings running with the field cutoff. 

4. Optimization and Interpolation 

At a given order K and coupling A, we can adjust 4> ma x{^,K) in order to 
minimize or eliminate the discrepancy with the (usually unknown) correct 
value. As (f> ma x is varied, the curve (or the derivative) of the approximate ex- 
pression crosses the numerical curve (or its derivative). The strong coupling 
can be used to calculate approximately this optimal ma a;(A, K) and so this 
is a natural approach for the interpolation between the weak and strong 
coupling regimes. This procedure has been illustrated with two examples 7,8 . 
The calculation of the modified coefficients at {4>max) fall in three categories. 
Low k (the usual ones with exponentially small corrections calculable semi- 
classically) ; intermediate k (crossover; complicated but with universal fea- 
tures) and large k (power suppressed; no k\ behavior). The intuitive flow 
picture is that the beginning of the series corresponds to the behavior of the 
scaling variables near the gaussian fixed point. The large order, corresponds 
to the approach of the high-temperature/strong-coupling fixed point. The 
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coefficients in the crossover {4> max dependent) correspond to the crossover 
in the flows. 

5. Lattice Gauge Theory with One Plaquette 

A simple example where the general ideas discussed above can be easily 
implemented is a SU(N) lattice gauge theory with one plaquette 8 . The 
partition function reads. 



lep 

with (3 = 2N/g 2 . Specializing to N = 2, fixing the gauge to the identity for 
three of the links and integrating over the angular variables, we obtain 



Note that that due to the compactness of the group, the partition function 
has a large field cutoff, and also a large action cutoff which is gauge invari- 
ant. To construct a weak coupling expansion, we expand the square root in 
power of /3 -1 . However, the "coefficients" depend on (3 because the range 
of integration does. To get a regular series, we add the tails of integration 
because they are e~ 2/3 = e~ 8 / g effects, but this affects the asymptotic be- 
havior of the series and the coefficients now grow factorially. The optimal 
order to truncate this series 2/3—1/2. Incidently, this is also the order where 
the peak of the integrand moves out of the range of integration if it is kept 
below 2/3 as in the exact expression Eq. (3). 

It is possible 8 to keep the finite bound of integration or to modify it 
in order to minimize the difference between the series at a given order and 
the original integral as explained in the previous section. However, if the 
order is large enough, the optimal cutoff becomes close to 2/3. Note also 
that if wc simply consider the regular perturbative series truncated using 
the rule of thumb, it is easy in this simple example to define and calculate 
the non-perturbative part of the integral. It consists in the higher order 
terms (to be calculated with the finite range of integration) minus the tails 
of integration that we have added. 

6. The Non-Perturbative Part of the Plaquette 

In the previous example, we can obtain a converging expansion by calculat- 
ing the coefficients of the /3 _1 expansion keeping a finite range of integra- 
tion. When there is more than one plaquette, this is much more complicated 






(3) 
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and we are presently developing practical methods to perform these calcu- 
lations. In the rest of this talk, we discuss the regular perturbative series of 
the average plaquettc of lattice QCD in the quenched approximation. We 
define the average plaquette 

P{(3) ee (l/Af p ) - (l/N)ReTr(U p ))j . (4) 

with J\f p the total number of plaquettes. The series of P has been 
calculated 9 ' 10 up to order 16. Despite the long series available, the factorial 
growth is not apparent. Instead, the series appears more like it has finite 
radius of convergence n > 12 . This is impossible because P takes different 
limits 13 when g 2 — > 0^. The only plausible explanation is that the partition 
function has a zero in the complex (3 plane near (3 ~ 5.7. Zeroes have been 
found 14 for imaginary values of (3 within predicted bounds 12 (and not 
found below the lower bound), however their numerical significance remains 
to be established . 

A simplified model 15-17 that is capable of producing a scries with co- 
efficients growing factorially is 

PaK die-?* (1 - t 33/167T 2 )- 1 - 204 / 121 . (5) 
Jti 

The new parameter (3 is related to the lattice (3 by a relation of the form 
/3 = /3(l + di //?+...). When deriving this expression as a sum of bubble 
diagrams, one realizes that t\ = corresponds to momenta at the UV cutoff 
and ti = 167r 2 /33 corresponds to the Landau pole. 

Note that the expression is quite similar to the one plaquette integral 
discussed in the previous section. This could in principle be compared with 
what would be obtained from the probability distribution for one plaquette 
after integrating over all the other links. Note also that in his ITEP lec- 
tures, M. Shifman emphasizes that it is necessary to introduce the gluon 
condensate in order to keep ti low enough and regularize the perturbative 
series, exactly as we advocated in lattice gauge theory. 

Expanding (1 — t 33/167r 2 )~ 1 ~ 204 / 121 in powers of t and extending the 
integration range to oo, one finds at leading order that the coefficients of 
the expansion grow like (33/167r 2 ) fe r[fc + 204/121]. According to the 
rule of thumb, we should truncate the series at an order 4.79/3 — 2.19. For 
(3 = 6 and d\ small this means an order of about 25. For d\ = —3, this order 
is lowered to 12. Using Eq. (1) and shifting to the lattice (3, we conclude 
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that the error is proportional to 

( -^204/121-l/2 c -(16 7 r 2 /33)/3 _ (g) 

Except for the power -1/2, this expression depends on (3 as the fourth power 
of the two-loop renormalization group invariant scale. In practice, the de- 
pendence on the power of (3 is quite weak in the region of /3 where an 
empirical envelope for the accuracy curves can be seen. In Fig. 2, the thick 
line has been drawn using the one-loop formula formula with an adjustable 
normalization constant 1.3 x 10 10 x e~( 16Tr / 33 ^ and looks like a decent enve- 
lope. When the order increases, will the accuracy curves reach the straight 



quenched QCD 
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Fig. 2. Difference between the perturbative series at order 1 to 10 and the numerical 
value of the plaquette (in Log 10 scale). As the order increases, the curves get darker. 
The thick curve is discussed in the text. 

line or will they maintain some curvature? Higher order extrapolations 11,10 
favor the second possibility with a power 4 of the force scale. Lower power 
of this scale provide good fits of the accuracy curves at various orders. This 
question will be discussed into more detail in a forthcoming publication 18 . 

7. Conclusions 

A better control of perturbative series of the standard model is necessary. A 
field cutoff drastically improves the asymptotic behavior of series. In simple 
examples, the field cutoff can be chosen to minimize the discrepancy with 
the uncut theory. There seems to be a connection between the crossover 
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behavior of the perturbative coefficients and crossover behavior of the RG 
flows. Numerical methods need to be developed to implement large field or 
large action cutoffs. Analytic methods remain to be developed to estimate 
the various parameters determined empirically while studying the difference 
between the perturbative series and numerical values in lattice QCD. 
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